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Abstract 

The consistent description of the nuclear response at low and high mo- 
mentum transfer requires a unified dynamical model, suitable to account for 
both short- and long-range correlation effects. We report the results of a 
study of the charged current weak response of symmetric nuclear matter, 
carried out using an effective interaction obtained from a realistic model of 
the nucleon-nucleon force within the formalism of correlated basis functions. 
Our approach allows for a clear identification of the kinematical regions in 
which different interaction effects dominate. 
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1. Introduction 

Short-range nucleon-nucleon correlations have long been recognized to 
play a prominent role in the electromagnetic response of nuclei. Comparison 
between the longitudinal responses at momentum transfer |q| > 1.5 fm~^, 
measured by inclusive electron scattering experiments, and the predictions 
of independent particle models shows that the data in the region of the 
quasi- free peak are sizably overestimated [l], 0]. On the other hand, the 
results of approaches based on Nuclear Many-Body Theory (NMBT) and 
realistic models of the NN interaction 0, 0], in which the effects of short- 
range correlations are taken into account, provide a satisfactory description 
of a large body of data j^. 

Long-range correlations, leading to the excitation of collective modes, 
are also known to be important. They become in fact dominant in the 
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region of low momentum transfer, in which the space resolution of the probe, 
A ~ |q|~^5 becomes much larger than the average NN separation. Most 
theoretical studies of this kinematical regime have been carried out within 
the Tamm Dancoff (TD) [61] and Random Phase Approximation (RPA) [3], 
using effective particle-hole interactions. 

The development of a theoretical scheme allowing for a consistent treat- 
ment of short- and long-range correlations requires the construction of effec- 
tive interactions based on realistic models of the NN force. This goal has been 
pursued by the authors of Ref. [81], who used the formalism of correlated basis 
functions (CBF) and the cluster expansion technique j^, to obtain an 
effective interaction from a truncated version of the state-of-the-art Argonne 
Vis potential The resulting effective interaction and the associated ef- 
fective transition operators have been applied in a systematic study of the 
weak response of nuclear matter j^. 

The development of a unified description of the nuclear response to weak 
interactions at low and high energy is important, its quantitative understand- 
ing being needed in a variety of different contexts. The region of neutrino 



energy E^, ~ 0.5 — 3 GeV [l^ is relevant to the analyses of long-baseline 



neutrino experiments, such as K2K and MiniBOONE (see, e.g., Ref. [13 
and references therein), while nuclear interactions of low energy neutrinos, 
carrying energies of the order of tens of MeV, are believed to determine 



supernovse evolution and neutron star cooling [15 



The effective interaction proposed in Ref.|8(] has been recently improved 



[mp 

with the inclusion of the effects of three- and many-nucleon forces 16|. The 
resulting interaction has been used to calculate the nucleon-nucleon scatter- 
ing cross section in nuclear matter, needed to obtain the transport coefficients 



from the Boltzmann-Landau equation |16l . 117 . 

In this paper, we report the results of a study of the charged current 
weak response of symmetric nuclear matter, carried out using the effective 



interaction of Ref. [16| . The main purpose of our work is identifying different 
interaction and correlation effects, all described within the same dynami- 
cal model, and determining the kinematical regions in which they play the 
dominant role. 

2. Formalism 

In the kinematical regime in which the non relativistic approximation is 
applicable, the nuclear response to a weak probe delivering momentum q and 
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energy uj can be written in the form 

^(^'^) = 1^E (0|Ot(q)|n)(n|0(q)|0)5(a; + Eo-i5;„) . (1) 

n 

where iV is the particle number and |0) and |n) denote the initial and final 
nuclear states, respectively. The transition operators are obtained expanding 
the weak nuclear current in powers of Iql/^ri, m being the nucleon mass. 
In the case of charged current interactions, at leading order one finds the 
Fermi and Gamow- Teller operators, whose expressions in coordinate space 
are {R = (ri, . . . ,rAr)) 

N N 

{R'\0^{q)\R) = 6{R-R') gvY, e^'^^'V =5{R-F^) Of(q) , (2) 

i=l 1=1 
N N 

{R'\0^^{q)\R) = S{R-R') qaY. e^''^»o-,r+ = S{R-R') ^ ©^(q) . (3) 

i=l i=l 

In the above equations, rj specifies the position of the i-th particle, cTj de- 
scribes its spin and is the isospin rising operator. 

In our formalism nuclear matter is described using correlated states, ob- 
tained from the Fermi (FG) gas states through the transformation j^, [lo| 

where \n) is a determinant of single particle states representing nonin- 
teracting nucleons at density p = 2kp/37['^, kp being the Fermi momentum. 
The operator F, embodying the correlation structure induced by the non 
perturbative components of the NN interaction, is written in the form 

N 

F{1,...,N)=S n (5) 

j>i=i 

where S is the symmetrization operator, accounting for the fact that, in 
general, [fijjik] 7^ 0. 

The two-body correlation function fij features an operatorial structure 
refiecting the complexity of the NN potential: 

fij = ^ ifTsirij) + 5slftT{r^j)S,j] Pts ■ (6) 

TS 
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In the above equation Vij = |rj — rj\, Pts is the operator projecting onto 
two-nucleon states of total spin and isospin S and T, respectively, and 
Sij = 3{(Ti ■ Yij){(Tj ■ Vij)/rfj — (cTj ■ (Tj). The shapes of the radial 
functions fTsi^ij) and ftrifij), are determined through functional minimiza- 
tion of the expectation value of the nuclear hamiltonian in the correlated 



ground state, carried out at the two-body level of the cluster expansion [18 

The CBF formalism naturally leads to the appearance of an effective in- 
teraction, acting on the FG basis states. At lowest order of CBF perturbation 
theory Ves is defined by 

where H = T + V is the nuclear hamiltonian, T is the kinetic energy operator 
and = ^^.^j fjj, tijj being the NN potential. 

As the above equation suggests, in principle the approach based on the 
effective interaction allows one to obtain any nuclear matter observables us- 
ing perturbation theory in the FG basis. However, the calculation of the 
expectation value of the hamiltonian in the correlated ground state, needed 
to extract Vcs from Eq.Q, involves severe difficulties. 

In this work we follow the procedure developed in Ref. [sl, whose authors 
derived the expectation value of Ves carrying out a cluster expansion of {H) 
of Eq.d?]), and keeping only the two-body cluster contribution. The resulting 
expression is 



(0|Fefr|0) = J](U|/12 



m 



\^3)a , (8) 



where the laplacian operates on the relative coordinate and the subscript a 
indicates that the two-nucleon state \ij)a is antisimmetryzed. 

We have computed the effective potential defined by the above equation 



using the truncated form of the Argonne f is potential referred to as w § [19 
Three-nucleon forces, whose inclusion is needed to reproduce the saturation 
properties of nuclear matter, have been described following the approach of 
Ref. 0, 

in which the main effect of three- and many-body interactions is 
taken into account through a density dependent modification of Vij at inter- 
mediate range. The Euler Lagrange equations for the correlation functions 



fxs and ftr have been solved using correlation ranges taken from Ref. [21 
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The calculation of transition matrix elements between correlated states 
is also non trivial. While in the FG model the Fermi and Gamow- Teller 
operators can only induce transitions to one particle-one hole (Iplh) states, 
in the presence of correlations more complex final states yield nonvanishing 
contributions to the response. In addition, the numerical determination of 
the matrix elements requires the use of the cluster expansion technique. 

In this work we only include transitions to correlated Iplh states and 
approximate the transition matrix element according to 

= (p^lfti-tK^^ftFlO)./'^ " (P^I(l+E3.)O(l+E3.,)|0) . (9) 

where Qij = fij — 1 and O = Yli Oi, Oi being the Fermi or Gamow- Teller 
operator (see Eqs.([2]) and ([3])). 

The authors of Ref . [8j used a different truncation scheme to obtain Mp^ 
from the cluster expansion formalism. Our choice is motivated by the require- 
ment of consistency with the definition of the effective interaction, Eq.(IH]). 
However, we have verified that the difference between the numerical results 
obtained from the two different prescriptions never exceeds few percent. 

3. Effects of short-range correlations 

The calculation of Mp^ has been performed on a cubic lattice, with a 
discrete set of Nh states specified by the hole momentum hi, satisfying the 
conditions hi < kp and hi -|- q > kp- The size of the basis has 
been determined requiring that the response of a system of noniteracting 
nucleon computed on the lattice agreed with the analytical result of the FG 
model. Obviously, the response obtained from the discrete set of final states 
consists of a collection of delta function peaks. A smooth function of a; 
has been obtained using a gaussian representation of the energy conserving 
(5-function of finite width a. For sufficiently small values of a, the results 
become independent of cr. 

Figure [T] shows a comparison between correlated and FG matrix elements 
for the case of Fermi transitions. Note that the effect of the sizable quenching 
produced by short-range correlations is enhanced in the calculation of the 
response, whose definition involves the correlated matrix element squared 
(see Eq.([T])). Similar results are obtained for the Gamow- Teller transitions. 

The calculation of the response requires the knowledge of the energies of 
the initial and final states. In the case of Iplh final states, the argument of 



5 



1.05 



1.00 



0.95 - 



0.90 



0.85 



0.80 




Figure 1: Fermi transition matrix element at |q| = 0.3 fm~^ (Eq.Q) as a function 
of the magnitude of hole momentum |h|. The calculation has been carried out 
using a basis of 3040 ph states. The dashed horizontal line corresponds to the 
result of the FG model. 



the energy conserving 5-function appearing in Eq.([T]), reduces to En — Eq = 
Cp — Ch, Ck being the energy of a nucleon of momentum k in nuclear matter. 
Within our approach, Ck can be calculated using the effective interaction and 
the Hartree-Fock (HF) approximation: 

ek = 7^ + Yl {hk\V,s\hk)a = ^ + Uhf{\^) ■ (10) 

\\i\<kp 

The numerical results turn out to be in close agrement with those of Ref. 
[i^ . obtained within the CBF approach using the Fermi-Hyper-Netted-Chain 
(FHNC) summation technique and a realistic nuclear hamiltonian. 

Figure [2] shows the nuclear matter response for the case of Fermi tran- 
sitions. Panels (A), (B) and (C) correspond to the kinematical regions 
|q| < kp^ kp < |q| < 2kp and |q| > 2kp, respectively. The solid lines 
represent the results of the FG model, while the squares show the response 
obtained within the correlated HF approximation, i.e. using the correlated 
matrix elements of Eq.Q and the energy spectrum of Eq. ffTOl) . In order to 
clearly identify the ~ 30 % quenching due to short-range correlations, whose 
size is consistent with the results of Fig. [1], we also show by diamonds the 
response obtained replacing the HF spectrum of Eq. ffTOj) with the kinetic 
energy spectrum, i.e. neglecting Uhf(M) in Eq. flTUl) . It appears that interac- 
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tion effects leading to the appearance of the mean field also play a significant 
role, pushing strength to energies well beyond the kinematical limit of the FG 
model. Panels (A), (B) and (C) correspond to momentum transfer |q| =0.3, 
1.8 and 3.0 fm~^ and basis size Nh =3040, 2637 and 2777, respectively. The 
same qualitative pattern is obtained in the case of Gamow- Teller transitions. 



4. Effects of long-range correlations 

Equation ^ shows that in the correlated HF approximation the Fermi 
and Gamow- Teller transition operators are replaced by effective operators, 
acting on FG states, defined through 

{ph\{l + Y,g^^)0{l + Y.9^M = {ph\0.m . (11) 

j>i j>i 

Note that the above effective operators are different from those of Ref. [8], as 
we have chosen a different truncation scheme of the cluster expansions. 
The FG ph states, while being eigenstates of the HF hamiltonian 

Hhf = 5Z ek , (12) 

k 

with Ck given by Eq. (fTOil . are not eigenstates of the full nuclear hamiltonian. 
As a consequence, there is a residual interaction V^es that can induce transi- 
tions between different ph states, as long as their total momentum, q, spin 
and isospin are conserved. 

We have included the effects of these transitions, using the Tamm Dancoff 
(TD) approximation, which amounts to expanding the final state in the basis 



of one Iplh states according to [23 



I/) = |q, TSM) = J2 c^''\P^h^, TSM) , (13) 

i 

where Pi = hj + q, S* and T denote the total spin and isospin of the particle- 
hole pair and M is the spin projection. 

At fixed q, the excitation energy of the state |/), cu/, as well as the 
coefficients cj^^'^ , are determined solving the eigenvalue equation 

H\f) = {HHF + V,,,)\f) = {Eo + uj)\f) , (14) 
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where Eq is the ground state energy. Within our approach this amounts to 
diagonahzing a x matrix whose elements are 



i/;^*^ = (Eo + ep, - ehJ5., + (Kpi, TSM\V,s\hjP„ TSM) , (15) 
In TD approximation, the response can be written as 

2 

5(0. -..r^o, (16) 



TSM n=l 



Nh 

J2{cl''%{h,p., T5M|Oeff(q)|0) 

1=1 



where {c^^^^)i denotes the i-th component of the eigenvector belonging to 
the eigenvalue ujn^^^ ■ 

The diagonalization has been performed using a basis of ~ 3000 ph 
states for each spin-isospin channel. The appearance of an eigenvalue lying 
outside the particle hole continuum, corresponding to a collective excitation 
reminiscent to the plasmon mode of the electron gas, is clearly visible in panel 
(A) of Fig. [3], showing the TD response at |q| = 0.3 fm~^ for the case Fermi 
transitions. For comparison, the result of the correlated HF approximation 
is also displayed. Note that the sharp peak arises from the contributions of 
particle- hole pairs with S = 1 T = 0. 

In order to identify the kinematical regime in which long range corre- 
lations are important, we have studied the TD response in the region 0.3 
< |q| < 3.0 fm~^. The results show that at |q| <1.2 fm~^ the peak corre- 
ponding to the collective mode in the S = 1 T = channel is still visible, 
althoug less prominent. However, it disappears if the exchange contribution 
to the matrix element of the effective interaction appearing in the rhs of 
Eq. ffTSl) is neglected. 

The transition to the regime in which short-range correlations dominate is 
illustrated in panels (B) and (C) of Fig. [31 showing the comparison between 
TD and HF responses at |q| = 1.5 and 2.4 fm~^, respectively. 

At |q| =1.5 fm~^ the peak no longer sticks out, but the effect of the mixing 
of ph states with S = 1 and T = is still detectable, resulting in a significant 
enhancement of the strength at large u. At |q| =2.4 fm~^ the role of long 
range correlations turns out to be negligible, and the TD and correlated HF 
responses come very close to one another. The calculation of the response 
associated with Gamow- Teller transitions shows a similar pattern. 
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5. Conclusions 



The CBF formalism employed in our work is ideally suited to construct 
an effective interaction starting from a realistic NN potential. The resulting 
effective interaction, which has been shown to provide a quite reasonable 



account of the equation of state of cold nuclear matter [16(], allows for a 



consistent description of the weak response in the regions of both low and 
high momentum transfer, where different interaction effects are important. 

The results of our calculations suggest that in addition to the HF mean 
field, which moves the kinematical limit of the transitions to Iplh states 
well beyond the FG value, correlation effects play a major role, and must be 
taken into account. While at |q| < 0.5 fm~^ long-range correlations, leading 
to the appearance of a collective mode outside the particle-hole continuum, 
dominate, at |q| 2.0 fm^^ the most prominent effect is the quenching due 
to short-range correlations. 

As a final remark, we emphasize that the results discussed in this paper, 
while being certainly interesting in their own right, should be seen as a step 
towards the development of a unified description of a variety of nuclear mat- 
ter properties, relevant to the understanding of neutron-star structure and 
dynamics 16|, UJ. ■ 
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Figure 2: Nuclear matter response for the case of Fermi transitions. Solid lines: FG 
model. Squares: correlated HF approximation. Diamonds: same as the squares, 
but with the HF spectrum of Eq. (llOp replaced by the kinetic energy spectrum. 
Panels (A), (B) and (C) correspond to |q| =0.3, 1.8 and 3.0 fm"-*^, respectively. 
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Figure 3: Nuclear matter response calculated within the TD (squares) and corre- 
lated HF (diamonds) approximations, for the case of Fermi transitions. Panels 
(A), (B) and (C) correspond to |q[ =0.3 , 1.5 and 2.4 fm~^, respectively. 
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